Magnetic Impurity Affected by Spin-Orbit Coupling: Behavior near a Topological 

Phase Transition 



o 

(N 
C 

00 



i 

CO 

c3 



i 

c 

o 
o 



o 
in 

o 



F. M. Hu, 1 T. 0. Wchling, 2 ' 3 J. E. Gubernatis, 4 Thomas Frauenheim, 3 R. M. Nieminen 1 
COMP '/Department of Applied Physics, Aalto University School of Science, 
P.O. Box 11100, FI-00076 Aalto, Espoo, Finland 
2 Institute of Theoretical Physics, University of Bremen, Otto-Hahn-Allee 1, D-28359 Bremen, Germany 
3 Bremen Center for Computational Materials Science, 
University of Bremen, Am Fallturm la, D-28359 Bremen, Germany 
4 Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

(Dated: January 29, 2013) 

We investigate the effect of spin-orbit coupling on the behavior of magnetic impurity at the edge 
of a zigzag graphene ribbon by means of quantum Monte Carlo simulations. A peculiar interplay 
of Kane-Mele type spin-orbit and impurity-host coupling is found to affect local properties such as 
the impurity magnetic moment and spectral densities. The special helical nature of the topological 
insulator on the edge is found to affect nonlocal quantities, such as the two-particle and spin-spin 
correlation functions linking electrons on the impurity with those in the conduction band. 



Introduction. Spin-orbit coupling (SOC) plays a cen- 
tral role in topological insulator (TI) materials it 
opens a gap inside the bulk but supports the gapless 
electron states on the boundary. On the boundary of 
TI, impurity potential scattering is restricted by time- 
reversal symmetry 0, Q : the backscattering among the 
electrons is allowed only if it is accompanied by a spin 
flip. This restriction gives rise to novel Kondo physics of 
magnetic impurities in TI 0, Q . Accordingly, the behav- 
ior of magnetic impurities in host materials with SOC 
has recently attracted interest. In two-dimensional (2D) 
systems, several theoretical papers 0, @] report that the 
presence of the SOC can in general protect impurity's 
magnetic moment from being totally screened. However, 
for the specific case of the Rashba coupling [H, [i| it has 
been reported that the SOC only makes a small or a 
high-order contribution to Kondo scaling. 

In general, SOC and energies associated with the 
Kondo screening of magnetic impurities can be of the 
same order of magnitude. It remains to be explored what 
effects arise when crossing over from weak to strong SOC. 
In particular, the question of how the impurity behaves 
when the host undergoes a topological phase transition 
from a normal state to TI is an open. 

In this paper, we consider an Anderson impurity on 
the edge of a zigzag graphene ribbon with a Kane-Mele 
type SOC gradually strengthened, so a topological phase 
transition occurs in the system. Using quantum Monte 
Carlo (QMC) simulations at finite temperature, we pre- 
form a comprehensive study on the properties of impu- 
rity. Calculating local physical quantities on the impu- 
rity site, such as the average double occupancy, mag- 
netic moment, spin susceptibility and spectral densities, 
we find that in general the SOC can support the local 
moment formation. We see that in the case of zero or 
weak SOC, impurity is mainly dominated by its coupling 
with localized edge states, realizing a situation similar 
to the zero band width Anderson (ZBWA) model. In 



the case of strong SOC, the edge states are delocalized 
and the impurity behaves as a spin in a normal metal. 
This difference between strong and weak coupling also 
manifests itself in the dependence of the local moment 
on chemical potential, which in graphene can be tuned 
by a gate voltage. We also study nonlocal linking of 
the electrons on impurity with those in the conduction 
band. In nonlocal two-particle correlations, we observe 
a set of distinct spin-momentum relations which show 
the helical locking and interplay between backscattering 
and spin flip. This interplay signals the appearance of a 
topological phase in the host. Additionally, we find the 
spin rotation symmetry in spin-spin correlation functions 
to be partially broken by the SOC. This finding agrees 
with a previous study of an Anderson impurity in 2D 
helical metal with variational method that is valid for 
large Coulomb interaction Q . In this work, we study the 
effect of electron-electron interaction over a wide region 
and document that the Coulomb interaction enhances the 
anisotropy in spin-spin correlation function and plays a 
complicated role in spin and momenta scattering. 

Model and methods. Our starting point is the Hamil- 
tonian H = H^-m + H± + Hi where -Hk-m is the Kane- 
Mele model for a zigzag edge graphene ribbon [ll|. It 
has two pieces: Hk-m = H t + H scn with H t being the 
usual nearest-neighbor hopping of tight-binding model 
in graphene 
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H so is the SOC term and s z is the z Pauli matrix, 
thus has opposite signs for opposite electron spins. 



The 



parameters v\. 



±1 depend on the orientation 



of the two nearest neighbor bonds as the electron hops 



2 



from site i to j: Vij = +1 if the electron makes a left 
turn to the second bond. It is negative if it makes a right 
turn. Hi is the impurity Hamiltonian 



Ud^d^d^d^ . 



Here Ed is the energy of the impurity orbital and U 
is Coulomb repulsion inhibiting the simultaneous occu- 
pancy of the orbital by two electrons. Finally H2 de- 
scribes the hybridization between impurity and one of 
atoms on the edge (located at A-sublattice site R a o), 



Our principal computational tools are the single- 
impurity QMC algorithm [l2| for computing the lo- 
cal thermodynamic properties of the impurity and the 
method of Bayesian statistical inference for comput- 
ing spectral densities. The QMC naturally returns 
the imaginary-time Green's function Gda{T > 0) = 
(d a (T)d\) of the impurity. With this Green's function, 
we can easily compute the magnetic quantities on the 
impurity site such as the expected values of magnet mo- 
ment square m 2 d = ((d\.d^ — djdj.) 2 }, the doubly occu- 
pancy ndfUdi — (d^d^-d^di) and the static impurity spin 
susceptibility 



dT(m d (T)m d (0)), 



(1) 



where (3 = T~ l , m^ir) — e TH rrid(0)e~ TH . Computing 
the imaginary-time Green's function also enables us to 
compute the spectral density A(uS) = ^2 a A a {uj) by nu- 
merically solving [lU 
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FIG. 1: (color online), (a) Fhe LDOS as a function of energy 
E/t with different values for A. From top to bottom A = 0, 
0.02*, 0.04*, 0.06*, 0.08*, 0.1* and 0.2*. (b)-(c): n dup n ddovn , 
m d , and x versus fi/t. Here, U = 1.2t, V = 0.65*, Ed = —U/2, 
and T _1 = 32t _1 . 



occupancy. We see that increasing the SOC results in 
the average double occupancy n^ridi decreasing because 
the SOC suppresses the effective hybridization between 
impurity and edge states. Near the Dirac point, the av- 
erage local magnetic moments (both m d and x) are en- 
hanced by the SOC. We also note that with weak SOC, 
i.e., 0.0 < A < 0.05*, the maxima of m 2 d and \ are n °t 
at fj, = while for the large A, they are at fx = 0. We 
comment that with small SOC the localized states on the 
edge are antiferromagnetically coupled to impurity states 
below the Fermi energy, so shifting fi from zero can de- 
couple these oppositely-aligned spin states and lead to 
the development of a magnet moment [l3| . Contrary to 
the weak SOC case, strong SOC delocalizes the states 
and induces a flat LDOS on the edge, making the LDOS 



Local properties. To gain insight about the conse- 
quences of the SOC, we solved the Kane-Mele model 
on zigzag graphene ribbon numerically and obtained the 
localized density of states (LDOS). For the LDOS in 
Fig.QTa), we see that when A = 0, a sharp peak exists at 
the Dirac point (E/t = 0) suggesting the existence of a 
strongly localized state on the zigzag edge. When A 7^ 0, 
the central peaks becomes smooth, and at A = 0.2*, their 
values near E/t = approach a constant and thus near 
the Dirac point become similar to the LDOS in a nor- 
mal metal. The SOC suppresses the LDOS and therefore 
weakens the localized states on the edge. 

In Fig.[TJ we show the QMC results of rid^ridi, m d and 
X versus the chemical potential /1 which can be tuned 
by an electric field. With Ed = —U/2 fixed and fi = 
the impurity site is half filled, and thus the magnetic 
moment is driven by avoiding the possibility of double 
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FIG. 2: (color online). The spin susceptibility x as a function 
of T for different values of SOC. U = 1.2*, V = 0.65* e d = 
— U/2 and /i = 0. The solid line presents the results for ZBWA 
model. 
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FIG. 3: (color online). The spectral density A(uu) versus u/t. 
(a) The SOC A = 0. The inset (al) shows the position (abso- 
lute value) of two inner peaks versus 1/U, and the inset (a2) 
shows the height of these two peaks versus 1/U 2 . (b) The 
SOC A = O.lt. In both figures, V = 0.65t, fi = 0, e d = -17/2 
and T _1 = 32*" 1 . 



near the Fermi energy similar to that of a normal metal. 
Thus we see that half filling (hole-particle symmetry) op- 
timizes the magnetic moment at /i = 0. 

Figure [5] shows the impurity spin susceptibility x as 
a function of temperature for different values of A. We 
fix hole-particle symmetry by setting n = and £d = 
—U/2. Doing so means the average electron occupancy 
of the impurity site is one; i.e., it is half filled. From the 
figure, we also see that with small < A < 0.05t and 
the lowing of temperature, the spin susceptibility first 
increases, then decreases and finally saturates (totally 
screened). But when 0.05i < A < O.lt, x nrs t increases 
and then directly becomes saturated. We note that when 
the SOC is gradually opened, x goes cross over form the 
behavior in ZBWA model (solid line in Figured]) to that 
in a normal metal [14[ . We propose this transition occurs 
because of the edge states delocalized by the SOC. 

In Fig. [3l we show additional spectral densities A(ui) 
for the impurity site. Here we have A(uS) = A{—uS) due to 
hole-particle symmetry. In Fig. [3^ arc results for A = 0. 
In this case the impurity states are strongly correlated 
with the edge states. There are four peaks. The inner 
two are from the localized edge state which is split by 
the impurity state. The outer two peaks are from the 
impurity level which is split by the Coulomb interaction 
into Ed and Sd + U. This four-peak structure is similar 
to that in ZBWA model. It is known that in ZBWA 
model the separations and weights of two inn er p eaks 
are proportional to U" 1 and U~ 2 , respectively [14J, but 
here we cannot see this dependence on U. 

In Fig. 13b, we show A(oj) for A = O.lt. Clearly vis- 
ible is the central Kondo resonance and the heights of 
two inner peaks decreasing with increasing A. This be- 



havior is consistent with the expectation that increasing 
A weakens the localized states. We also see that when 
we increase U, the height of the central Kondo peak, as 
well as that of the two inner peaks, decrease, but the 
heights of the outer two peaks from impurity levels in- 
crease. These behaviors are similar to those of an impu- 
rity in normal metal. Additionally, the two smooth peaks 
around u) = ±t are remnants of van Hove singularities. 

Nonlocal properties. The most interesting property 
of a TI boundary is time-reversal invariance caus- 
ing backscattering always being accompanied by spin- 
flip scattering. In order to directly see signatures 
of this phenomenon, we computed a set Pk of two- 
particle correlation functions between the impurity 
site and the conduction electron 15 with momentum 
k: P k = {(ck,td\.d±&_ k ^), (c-k^d^dicl^), (cfc^dj^cj^), 

{c-k,fdtdic_ k I )}. In Fig. [4j we focus on the states 
with k = ±0.001257r that are below Fermi energy. In 
Fig. UK, we see that for positive A, the correlation func- 
tion (cfe^djd^clj. j_) = while (c-k,f d|dj,c|. j_) ^ 0, and 
with negative A, the situation is opposite. This means 
that the conduction electrons with a fixed spin orienta- 
tion can only move in one direction. When the sign of 
SOC changes, they then move in the opposite direction. 
Backscattering occurs if accompanied by a spin flip. In 
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FIG. 5: (color online) 

(cfc,t4>( d J- c -fc,4.) and B1 is ( c -fc.t4)^ c fc,4.) 

(cfe,t«4)< d 4-4,4.) and D1 is ( c -k,t d l)( d l c -k,i) 
ters are the same to those in Fig. [4] 
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FIG. 6: (color online). The ratio Sf/Sf as a function of 
the position i on A-sublattice edge, (a) The non-doping case 
/i = 0, (b) the doping case /j. = — O.lt, and (c) spin-orbit 
coupling A = 0.15t and the chemical potential is varied. In 
(a)-(c) U = 1.2t; in (d) fj, = -0.2t, X = 0.154 and U is varied. 
In all figures V = 0.65t, e d = -U/2 and T" 1 = 32t _1 . 



Fig. 2b, we observe that the two correlation functions 
(ck^d^dicl j_) and (c-k,fd^.d±c_ k ^} vanish, meaning that 
spin flip combined with forward scattering is forbidden 
in helical liquid phase. 

To evaluate the effect of electron-electron correlation, 
we also computed the set Pk of products of two single- 
particle correlators, each corresponding a two-particle 
correlator P k ; for example, we computed the product 
(ck,fd^) (dic_ k ^) that is part of (ck,\dLd\,c*_ k ^). The dif- 
ference between the two corresponding functions is just 
the vertex contributions of interaction [17J . We show the 
products of the two single-particle correlation functions 
in Fig. [5j and comparing Figs. |4] and we can see the 
helical contributions to the two-particle correlations. 

If we study the effect of the Coulomb interaction by 
varying U over a wide range and fixing Sd = —U/2, the 
difference between a member of Pk and the correspond- 
ing member of Pk becomes more obvious. Figure [5] shows 
that the product of the two single-particle correlators 
monotonically decreases as U increases, while the two- 
particle correlators at first increase (0 < U < 1.2t) and 
then decrease (U > 1.2t). In fact we can enhance i\ by 
increasing e^ — Ek}^ 1 (Ek is the energy of ±fc states) and 
the impurity's double occupancy. U decreases these fac- 
tors. While a member of Pk not only depends on coupling 
between the impurity level Sd and ±/c states but also de- 
pends on the impurity's moment. U at first develops the 
magnetic moment and then suppresses the SOC. 

Using the same algorithm as those used for computing 
Pk and Pk, we computed the spatial spin-spin correla- 
tion functions, Sf, Sf and Sf defined as Sf = {(dld^ — 

4^)( c iW~ c l c u))> s f = ((4^+4 d t)(4t c u+ c ii c it)>> 



FIG. 7: (color online). The ratio Sf/Sf. Here i is a B- 
sublattice position on the edge. All the parameters are the 
same as those in Fig. [6] 



and Sf - H) 2 <(4<^ ~ 4 d t)(4r c U ~ c l c n))- We find 
that the spin rotational symmetry in these correlations 
is partly broken by the SOC, i.e., S? ^ Sf but S? = Sf, 
and that this asymmetry agrees with what was found for 
an impurity in the two-dimensional helical metal 0] . 

In Figs. [5] and [3 we display this anisotropy by plot- 
ting the ratio Sf / Sf on the two sublattices. It is clear 
that the behaviors on two sub-lattices are similar but in- 
dependent: the ratio equals one when A = 0, but when 
A ^ 0, the ratio Sf /Sf has a large distance asymptotic 
behavior ~ a sin(2/3i?j). Here Ri denotes a site on the 
zigzag edge. The two parameters a and (3 depend on the 
SOC, Coulomb interaction and chemical potential. The 
magnitude of a increases as A and U increase, meaning 
that SOC and electron-electron interaction enhances the 
anisotropy in spatial correlation functions. Further, (3 
equals 7T when /.t = or A = 0, and it decreases as A 
increases but increases as \x moves below zero. 

Discussion and conclusion. In this paper, we study an 
Anderson impurity on a zigzag graphene ribbon undergo- 
ing a topological phase transition driven by a Kane-Mele 
type SOC. Using QMC simulations, we investigate both 
the local and non-local properties of the impurity. We 
find that with the SOC increasing from zero, the for- 
mation of local moment is supported and the impurity 
behavior crosses over from a spin in a band with zero 
width to one in a normal metal. This is because SOC 
delocalizes states on the edge. As for non-local proper- 
ties, we clearly see the helical locking and the collabo- 
ration of backscattering and spin flip in the two-particle 
correlators and find a broken spin rotation symmetry in 
spin-spin correlation functions. 

Although the intrinsic SOC of carbon atoms is weak 
[3], a strong SOC can be achieved by heavy-atom 
doping 3- 23 1 and tuned by a gate bias 24 1 . The mag- 
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netic impurity could be realized by dangling a bonds in 
graphene [25|. Several physical quantities discussed can 
be experimentally measured: scanning tunneling micro- 
scope (STM) can measure A(cj), and recent developments 
in the field of spin-polarized STM [26| open the possibil- 
ity to detect the spin-spin correlations. 
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